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A SIMPLE APPROXIMATE METHOD FOR OBTAINING- SPANWISE 
LIFT DISTRIBUTIONS OVER SWEPT WINGS 
By Franklin W. DIederich 


SUMMARY 


It is shown, how Schrenk's empirical method of estimating the lift 
distribution over straight wings can he adapted to swept wings hy 
replacing the elliptical distribution hy a new "ideal" distribution 
which varies with sweep. The application of the method is discussed 
in detail and several comparisons are made to show the agreement of the 
proposed method with more rigorous ones. It is shown how first-order 
compressibility corrections applicable to subcrltical speeds may be 
included in this method. 


INTRODUCTION 


A great number of methods for estimating the lift distribution 
over unswept wings have been available ‘for some time, most of them based 
on the concept of the lifting line. The accuracy of these rational 
methods is limited by the degree of closeness with which the lifting 
line represents the physical conditions, so that much of the time 
consumed in obtaining veiy accurate mathematical solutions for the 
lifting-line equation Is unwarranted. However, idle experience gained 
with solutions obtained by these more accurate methods has pointed 
the way to simplified empirical methods which do not involve the 
solution of any mathematical equations. 

Glauert’s suggestion (reference l) that the shape of the lift 
distribution was "intermediate between that of the aerofoil (wing) and 
that of the ellipse" was followed up by Schrank (reference 2). He 
decided that in the case of untwisted wings the lift distribution 
could be represented approximately by the arithmetic mean of an ideal 
distribution and the plan-form distribution of equal area, the Ideal 
distribution being of elliptical shape. For twisted wings he resolved 
the lift distribution into an additional distribution to be obtained 
as described above and a basic distribution to be obtained by taking 
one -half of the distribution of the product of chord and angle of 
twist (measured from the angle of zero lift of the wing) and rounding 
off any sharp comers in that distribution. By means of several 
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numerical comparisons he demonstrated that the results of his method 
checked those of the rational methods and experimental results with 
sufficient accuracy for many practical purposes. Further comparisons 
made in references 3 and 4 corroborate this observation. 

The classical lifting-line methods do not apply to swept wings, 
however, so that Schrehk's approximate method, which is "based on the 
results of lifting-line methods, must be expected to fail, as well. A 
number of rational lifting-surface methods have been developed to treat 
swept wings (for example , references 5 and. 6 ), which are considerably 
more time-consuming than even the lifting-line methods. These lifting- 
surface methods have been used for a number of computations, and 
corroborative experiments have been carried out on a set of full-scale 
models (reference 7) • The results provide a basis for a method similar 
to Schrenk's, which is outlined in the present report. The only differ- 
ence consists of the fact that the ellipse is no longer considered as 
the ideal distribution for the additional lift distribution. A set of 
new ideal distributions which depend on sweep has "been determined as the 
set of curves which, when averaged with the chord distribution, yields 
results which are closest to theoretical lifting-surface solutions for 
various plan forms and degrees of sweep. The theoretical solutions 
required were available in a number of published and unpublished reports 
for angles of sweep up to ^ 5 °; no ideal distributions have been derived 
for sweep angles larger than 45° « The distributions are independent of 
the aspect ratio within the range of practical aspect ratios, as in the 
case of the unswept wing. They can be used for estimating the additional 
loading (and, consequently, the total loading) over any swept wing with 
a degree of accuracy which is entirely adequate for many purposes. 


SYMBOLS 


A aspect ratio (b 2 /S) 

b span, feet 

c chord, measured parallel to the plane of symmetry, feet 
c average chord, feet (s/b) 

Cl wing lift coefficient 

c 1 local lift coefficient 

f ordinate of the "ideal" distribution curve 

G Glauert-Prandtl correction ^l//l - M 2 ^ 

k5 angular change in zero-lift direction of any section produced by 
flap displacement 6 , radians 
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Mach n umb er 
slope of wing lift coefficient ((LCjJdjxJ 
slope of section lift coefficient 


(KMa=») 


average slope of section lift coefficient 
wing area, square feet 
lateral ordinate, feet 

local angle of attack measured from a common reference, radians 

average local angle of attack measured from the common reference 
(assumed to represent the zero-l if t angle of the wing), radians 

local angle of attack measured from the zero-lift angle of the wing, 
radians 

angle of sweephack at the quarter- chord line (negative values of A 
indicate sweepf orward) , degrees 

taper ratio (c tlp /c root ) 


Subscripts: 

additional lift distribution (for = l) 
b basic lift distribution ("for = 0 ) 
c compressible flow 

e equivalent (in incompressible flow) 

s equivalent unswept wing 

DE3CBZPTI0M OF TEE METHOD 
Additional Lift Distribution 

The additional lift distribution is obtained from the relation 


cc 


l 


cCp, 


1 


c 



( 1 ) 
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where f, which replaces the elliptic distribution — i/l - of 

* V kw?) 

references 2 , 3 , and 4, is to be read from figure 1 for the given sweep 
angle A, measured to the quarter-chord line. The value of Iq ia 
obtained from the definition 



lb/2 


®oC dy 


( 2 ) 


by graphical or numerical means, as described in references 3 a 111 ! 4. 

Wings of variable sweep, such as the M, W, and pterodactyl typeB, 
have been proposed from time to time. It appears that the lift distribu- 
tion over such a wing cannot be estimated very accurately by any simple 
means. One way of obtaining a useful approximation, however, appears to 
be the use of_ an ideal distribution f corresponding to the average 
sweep angle A defined by 


flb/2 

tan A = — tan Ac dy 

S. 

Uo 


(3) 


Basic Lift Distribution 


The basic lift distribution is less affected by sweep than the 
additional lift distribution. Furthermore, it usually constitutes a 
small part of the total lift distribution. Consequently, no essential 
deviation from Schrenl’s method appears to be warranted. The basic lift 
distribution is obtained from the relation 



™o c m 
2 5 ^ 0 


where 


and 


a Q = a - a 




ac dy 


( 4 ) 

(4a) 


(4b) 
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All s har p comers are rounded off in accordance with the suggestions of 
references 3 and k- in such a manner as to yield a total lift equal to 
zero. The factor m/®o ^ Las ^ eerL introduced "because comparisons of "basic 
lift distributions estimated by Schrenk ' s method with those calculated Ty 
lifting- surface methods indicate that much better agreement between the 
two can be obtained if Schrenk 's results are multiplied by the factor m/S^ 
this applies to both swept and unswept wings. 


Subcritical Compressibility Corrections 

Analyses of linearized three-dimensional compressible flows (for 
example, reference 8 ) indicate that the lift of a wing in a compressible 
flow can be obtained by finding the lift of an equivalent wing in 
incompressible flow and multiplying it by the Glauert-Prandtl correction G-. 
The aspect ratio and the sweep of the equivalent wing are related to these 
of the actual wing as follows: 



(5) 


tan Aq = G- tan A 


(5a) 


The taper ratio is the same for both wings. 

This procedure results in lift distributions over unswept wings 
which are unaffected hy compressibility. The lift distributions over 
swept wings, however, are modified by coup ress Ability In the. same 
manner as they are hy increases in the angle of sweep in incompressible 
flow. Specifically, the ideal distribution f for a wing in com- 
pressible flow is to be chosen for the angle A e rather than A. 

Similarly, the lift slope of a wing in compressible flow may be 
obtained from that of an equivalent wing in incompressible flow. An 
unpublished analysis which takes account of the effect of sweep on 
dovnwash in incompressible flow in a rational but approximate manner 
has "been performed by Mr. Thomas A. Toll. According to this analysis, 
the lift slope of a swept wing in incompressible flow may be estimated 
from the relation 


(A + 2) cosA 

m = — nij, 

A + 2 cos A 


(6) 


where the lift slope nig is to be taken for an unswept wing of the 
same aspect ratio and taper ratio. If the effects of compressibility 
are accounted for in the manner outlined previously, equation ( 6 ) is 
modified as follows: 

= ^ + g ) cos Aq _ 

Aq + 2 cos c (7) 


m 


y* 
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The lif t slope nig Is to "be taken for an im swept vlng of the actual 

c 

aspect ratio and taper ratio operating at the actual free-stream Mach 
number. It is approximately equal to G times the slope nig of a 
wing with the actual taper ratio and the aspect ratio Aq. 

These corrections are "based on the linearized equations of flow 
and consequently "become invalid aB the free-stream velocity approaches 
either that of sound or that corresponding to the critical Mach number 
of the wing. A more detailed a n al y sis must "be made if either a Mach 
number of about 0.9 or the critical Mach number of the wing (whichever 
is lower) is exceeded. 


DISCUSSION 


In considering the reliability of the method outlined in this 
report, it must be kept in mind that the method 1 b based on results 
obtained by means of potential-flow lifting-surface, methods and 
consequently does not take account of boundary-layer effects. Both 
this method and the more refined methods are useful because they 
provide first-order estimates of the load distribution at low and 
and medium angles of attack; they also provide bases for ' empirical correc- 
tions for boundary- layer effects where these effects are of particular 
importance, that is, for large angles of attack and sweep. Since, 
however, the boundary-layer effects depend on Beynolds number, Mach 
number, the section properties, and the plan-form parameters (which 
include sweep), the experimental information now available is entirely 
insufficient for the determination of a set of boundary-layer correc- 
tions. When more experimental information is available and such a set 
of corrections is obtained, it will apply to the approximate methods 
as well as to the rational potential-flow methods of estimating the 
load distribution. 

Figure 2 shows the agreement of the additional, lift distributions 
obtained for different wing configurations by the empirical method of 
this report with those calculated by Falkner’s method (reference 5 ) 
and thoBe measured on full-scale models (reference j). It appears that, 
while the results estimated by the two methods differ slightly from 
each other, both agree equally well with the experimental results, at 
least at the relatively low angles of attack used in the tests. 

Figure 2(d), which shows the agreement of the approximate method with 
Falkner’s in the cane of an untapered sweptback wing, indicates that 
the approximation Is valid even for large taper ratios. No experi- 
mental results appear to be available for comparison. 

Figure 2(e) has been drawn for a pterodactyl wing in an attempt 
to show how accurate a result may be obtained for wings of variable 
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sweep by tiLe approximate method. The Ideal distribution, f for the 
pterodactyl wing shown in the figure was chosen for an average sweep 
angle A = 25 0 . The agreement of the distribution obtained in this 
case with that calculated by Falkner’s method is adequate for nest 
structural purposes. 

The agreement between the results obtained by Schrenk’s method, 
the lifting-line theory (reference 9 ), and. Weissinger’s method (refer- 
ence 6 ) for flapped and twisted untapered wings of aspect ratio 5 is 
shown in figures 3 (a) and 3(b), respectively. In fairing the distribu- 
tions for Schrent’s method equal areas have been subtracted at the tip 
and at or near the root. It appears that the results of Schrenk’s 
method are slightly lower than those of lifting-line theory and that 
they would have to be reduced about 30 percent in order to agree with 
those of Wels 3 inger's method. Application of the factor m/mi^ would 
reduce the value by 32 percent; this fact was the basis for including 
the factor in equation (4). 

It appears that there Is comparatively little difference between 
the distributions for the three different values of sweep, although the 
usual trends are noted, in that the tendency of both sweepback and 
sweepf orward is to reduce the ordinates of the lift-distribution curve 
for a given angle of attack and the tendency of sweepback Is to shift 
the load outboard, whereas that of sweepforward Is to shift the load 
inboard, Inasmuch as these differences are smaller than those between 
the estimated distribution and the calculated distribution for zero 
sweep, and inasmuch as the basic lift distribution usually forms only 
a small part of the total distribution. It appears that little gain 
could be had by attempting to obtain a closer estimate for it. 

It miight be expected that of the family of curves presented in 
figure 1 the one for zero sweep would be an ellipse. Actu ally , there 
is a slight deviation, since the family of curves is based on lifting- 
surface theory, whereas Schrenk’s adoption of the ellipse is based on 
lifting-line theory. Experimental results for the total lift of wings 
with zero sweep agree more closely with those of lifting-surface theory 
than those of lifting- line theory (reference 6 ), and, while few corre- 
sponding comparisons appear to have been made for the lift distribu- 
tions on wings of zero sweep, the sane relative accuracy may be expected. 
A lifting-line calculation has been made for the wing of figure 2(a) 
by Multhopp’s lifting-line method (reference 9) and It appears that 
FaXkner’s distribution Is in slightly better agreement with experimental 
values than the lifting-line distribution. 

The fact that the ideal distribution curves do not pass through 
zero at the wing tip results from the consideration that in order to 
yield zero lift at the tip the positive ■ value of the tip chord requires 
an equal negative value of the ideal distribution to be averaged with 
It. Since the tip chords vary for different plan forms, it has been 
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found to "be expedient to draw the distributions only up to the point 
where they reach zero. If in using them any sharp comers are obtained 
in the lift distribution, they may be faired by eye. The areas under 
the ideal distribution curves have been so adjusted as to give an 

CC y y 

area of 1.000 under the curves for — _ plotted against — with an 

5Cl Vs 

allowance made for the fairing that may normally be required. 

The term "ideal" applied to the distributions of figure 1 has been 
carried over from reference 2 and is actually somewhat of a misnomer 
in the case of swept wings since, by Munk's stagger theorem, the 
elliptic distribution always causes the least induced drag, regardless 
of sweep. The term is not intended to imply that the distributions of 
figure 1 constitute desirable lift distributions. 


CONCLUDING EEMAEKS 


It has been demonstrated that the lift distribution over swept 
wings can be estimated with adequate accuracy for most practical 
purposes by means of simple extensions to Schrenk’s method. These 
extensions consist primarily of using a new set of ideal lift- 
distribution curves which depend on sweep instead of the ellipse in 
estimating the additional lift distribution. The shape of the basic lift 
distribution curves can be estimated by Schrenk's method without any 
modification for all practical sweep angles, but the magnitude should 
be multiplied by the ratio of the wing lift slope to the average section 
lift slope for all cases. Including the straight wing. 

The results of this empirical method are in as good agreement 
with experimental results as those furnished by the rational lifting- 
surface methods for all, practical sweep angles and for tapered wings. ■ 
For wings with sweep angles which vary along the span the agreement 
between the empirical and rational results is not as good, but may be 
adequate for structural purposes. 

Insufficient experimental Information exists at this time to 
provide a basis for empirical corrections to account for boundary- 
layer effects; this shortcoming affects the results of both the 
empirical and the rational methods. However, at low and medium 
angles of attack no correction Is required for most purposes. 

In order to estimate the shape of the distribution curve for 
compressible flows the Ideal-distribution curve must be chosen for 
an increased sweep angle. A first-order estimate of the magnitude 
of the curve or, for that matter, the lift slope is given by equation ( 7 ) 
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Both, these corrections are applicable only to subsonic and subcritical 
speeds. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 


RJiS’ERENCES 


1. Glauert, H. : The Elements of Aerofoil and Airscrew Theory. 

Cambridge Uhiv. Press, 1926* p. 155* 

2. Schrenk, 0. : A Single Approximation Method for Obtaining the 

Spanwise Lift Distribution. NACA TM No. 948, 1940. 

3 . Anon.: Airplane Airworthiness. A Single Approximate Method of 

Obtaining the Spanwise Distribution of Lift on Wings. Civil 
Aero. Manual 04, CAA, U.S. Dept. Commerce, Feb. 1, 194-1, 
appendix V. 

4. Flatt, J. : Evaluation of Methods for Determining the Spanwise 

Lift Distribution. ACTE No. 4952, Materiel Command, Army Air 
Forces, June 23 , 1943- 

5 . FaLkner, V. M. : The Calculation of Aerodynamic Loading on Surfaces 

of Any Shape. E. &. M. No. 1910, British A.B.C., 1943. 

6 . Weissinger, J.: The Lift Distribution of Swept -Back Wings. 

NACA TM No. 1120, 194?. 

7 . Van Dorn, Nicholas H., and DeYoung, John: A Comparison of Three 

Theoretical Methods of Calculating Span Load Distribution on 
Swept Wings. NACA TN No. 1476, 1947. 

8 . G 8 thert, B.: Plane and Three-DLmensional Flow at High Subsonic 

Speeds. NACA TM No. 1105, 1946. 

. Multhopp, H. : Die Berechnung der Auftrleb3verteilung von 

Tragfliigeln. Luf tfahrtf orschung, Bd. 15, Lfg. 4, April 6 , 

1938 , pp. 153-1&9* (Available as E.T.P. Translation No. 2392, 
British Ministry of Aircraft Production. ) 


9 



10 


NA.CA EM TTo. L7I07 







Lateral ordinate, ~- 
! % 


wnwL «M*orr 
««n« ran notowumoi 


Figure z.- Comparison of additional hit distributions for swept wings. 
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(a) Wing 'with 55.6 span flap . 







